Complementarity of Weak Lensing and Peculiar Velocity Measurements 

in Testing General Relativity 
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We explore the complementarity of weak lensing and galaxy peculiar velocity measurements to 
better constrain modifications to General Relativity. We find no evidence for deviations from GR 
on cosmological scales from a combination of peculiar velocity measurements (for Luminous Red 
Galaxies in the Sloan Digital Sky Survey) with weak lensing measurements (from the CFHT Legacy 
Survey). We provide a Fisher error forecast for a Euclid- like space-based survey including both 
lensing and peculiar velocity measurements, and show that the expected constraints on modified 
gravity will be at least an order of magnitude better than with present data, i.e. we will obtain ~ 5% 
errors on the modified gravity parametrization described here. We also present a model-independent 
method for constraining modified gravity parameters using tomographic peculiar velocity informa- 
tion, and apply this methodology to the present dataset. 
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I. INTRODUCTION 

The discovery of the cosmic acceleration [1, 2] a decade 
ago has forced cosmologists to modify their simple pic- 
ture of the Universe, i.e. a universe dominated only by 
matter and described solely by General Relativity (GR). 
Over the last ten years, physicists have suggested two 
main avenues for explaining the late-time acceleration 
of the Universe. The first is the introduction of Dark 
Energy (DE), with an effective negative pressure, which 
dominates the late-time energy density of the Universe 
thus causing the acceleration. An alternative explana- 
tion is to modify the law of gravity on large scales thus 
altering the predicted expansion history of the Universe 
to be in line with the observations. In this paper, we fo- 
cus on this latter explanation for the cosmic acceleration 
and present combinations of observables that can be used 
to test the validity of GR on cosmological scales. 

General Relativity is a metric theory of gravity that 
can describe the relationship between matter perturba- 
tions, gravitational potential and space curvature pertur- 
bations. Several authors have shown that by combining 
various probes of the large-scale structure in the Uni- 
verse, it is possible to test the relationship between these 
quantities which, in the linear regime, can generally be 
described by two functions of time and scale [3-10] (see 
[11] for a review and references therein). For example, 
several recent attempts have been made to constrain pa- 
rameters associated with these two functions using the 
latest weak lensing measurements [12-15] and these cur- 
rent analyses appear to be consistent with GR once sys- 
tematic errors are taken into account. However, these 
current constraints on Modified Gravity (MG) parame- 
ters are still weak, especially as there are degeneracies 
between the parameters that cannot be broken by weak 
lensing measurements alone. 

We show in this paper that these degeneracies can be 



broken through a combination of weak lensing (WL) and 
peculiar velocity (PV) measurements (see [16] for an ear- 
lier approach). This is motivated by the fact that the 
WL experiments probe the lensing potential, which de- 
termines the trajectories of photons through the Uni- 
verse, while PV measurements probe the gravitational 
potential that governs the dynamics of galaxies indepen- 
dent of their galaxy bias. A combination of these two 
observables can, in principle, simultaneously measure the 
MG parameters used to describe the possible modifica- 
tions to the relationships between the metric perturba- 
tions defined in GR [3, 4]. 

This paper is organized as follows. In Section II, we 
introduce the two MG parameters that describe the re- 
lationships between the GR metric perturbations. In 
Section III, we use current observational data to present 
constraints on these MG parameters. In Section IV, we 
perform a Fisher matrix error forecast for a future Euclid- 
like space-based DE mission, while in section V we dis- 
cuss a new method for extracting model-independent in- 
formation from PV measurements. We conclude in Sec- 
tion VI. 



II. MODIFIED GRAVITY PARAMETERS 

Linear metric perturbations around a background 
Friedmann universe are described by the line element in 
the Newtonian gauge 



-{l + 2^)dt'^ + {l + 2^)a' 



5ij dx^ dx-' 



(1) 



where $ and ^' denote the space curvature perturbation 
and the gravitational potential respectively. The dynam- 
ics of pressureless matter fluctuations are determined by 
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energy-momentum conservation, and are given by 

S+- = 0, (2) 
a 

e + HO = — (3) 
a 

where k is the wavcnumbcr of the perturbations, 5 is the 
energy density contrast, 6 is the divergence of the veloc- 
ity perturbations, and the overdot denotes the derivative 
with respect to the physical time t. In this paper, we 
only consider perturbations on sub-horizon, but linear, 
scales (i.e. H <$:k/a< 0.1 h/Mpc). 

Two additional equations, which need to be provided 
by a theory of gravity, are required to solve the evolution 
equations for the four scalar perturbations given above 
(<f>, 5*, 6 and 9). In this paper, we use the following modi- 
fied gravitational equations to parametrize possible mod- 
ifications from the expected GR relations between these 
perturbations (see [8], and references therein, for a dis- 
cussion of the various possible parametrizations) , 

fc2^ = ~47rGa^n{k,a)pS, (4) 
fc2($_^) = 8TrGa'^T.{k,a)pS, (5) 

where G is the Newton constant measured in a 
Cavendish-type experiment. Here, the function fi charac- 
terizes a modification of Newton's constant in both space 
and time, and the S function describes a modification to 
the lensing potential $ — again with space and time. 
Since /i affects the gravitational potential ^f, it changes 
the growth rate of density fiuctuations S via Eqns. (2) 
and (3). Weak lensing measurements probe the lensing 
potential, and are thus affected by both S and ^ (via S). 
Therefore there is a degeneracy in the constraints on S 
and /i from WL. 

On the other hand, peculiar velocities of galaxies are 
determined by the gravitational potential through the 
Euler equation, Eqn. (3), and are thus affected by fi but 
not by E. Therefore, the combination of WL and PV 
measurements can break the WL degeneracy between S 
and fj., and allow cosmologists to probe these two func- 
tions separately. Furthermore, [17] recently studied the 
theoretical priors on /i and E from scalar-tensor gravity 
theories, clustering dark energy models and interacting 
dark energy models, and showed that each of these theo- 
ries has a distinct path in the /i-S parameter space thus 
providing the opportunity to distinguish between these 
possible theoretical models. 

For the purpose of highlighting the key features of the 
combined constraints, in the following sections we use a 
simplified parameterization of S and ^ given by 



(6) 



where s specifies the power of the assumed time variation 
of S and ^, and Es and Hs are constants. In Section IV 
we will consider two cases, s = 1 (the linear model) and 
s = 3 (the cubic model). The linear model is motivated 



by DGP [18-20], while the cubic model is motivated by 
general plausibility arguments that fi may change in pro- 
portion to matter density [12, 21]. 

The parametrization in Eqn. (6) assumes that E and fi 
are scale independent. It should be emphasized that this 
is simply a restriction that arises when using the current 
data sets. It has been shown that the WL data is far 
more sensitive to scale-dependent modifications to GR 
[7, 8, 14]. Moreover, the PV observations we use in Sec- 
tion III assume that the growth rate is scale-independent 
[22], and therefore this data cannot be used for tests of 
scale-dependent models. In order to make a consistent 
comparison between current errors and future forecasts, 
we also do not consider the scale-dependent models in 
our future forecasts in Section IV. A comprehensive fore- 
cast of general scale and time dependent MG constraints 
from future PV measurements will be presented sepa- 
rately [23]. 



III. CURRENT CONSTRAINTS 

In this section, we present the current constraints on 
fi and S defined in Section II from the weak lensing 
measurements of the CFHT Legacy Survey (CFHTLS) 
[24, 25], and PV measurements obtained from Sloan Dig- 
ital Sky Survey (SDSS) Luminous Red Galaxies (LRGs) 
selected from their Data Release Seven (DR7) [22]. 



A. Weak lensing measurements 

The WL shear power spectrum is equal to the conver- 
gence power spectrum given by [6] 



367r 
'25 



dlnfcA: 



n 



dzW{z)n[kx{z)]{^-'f){k,z) 



where 
spectrum. 



A2 



(7) 



denotes the primordial curvature power 
x(z) is the comoving distance and the 
window function for WL is defined as W{z) = 
dz'n{z')[x{z') /x{z) — 1]. Given C^, we can obtain 
the E-mode component of the shear correlations via 
[26], 



(8) 



where Jq is the Bessel function of the first kind. 

For the WL observable, we use the data from the 
CFHTLS- Wide survey, presented in [24]; this is derived 
from shear measurements for 2 x 10^ galaxies with mag- 
nitudes 21.5 < Iab < 24.5 in a 35 sq. deg. effective area. 
We note that there is a modest systematic effect in the 
CFHTLS data, due to residual field-to-field variations in 
shear estimation on the scale of the camera field of view 
[14]. 

We estimate the data covariance matrix using the Hori- 
zon simulation as described in [27], section 3; we measure 
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FIG. 1: The constraints on Es and fj,s from the latest observational data; s — 1 linear model on the left and the s — 3 cubic 
model on the right. The inner dark-shaded contours are the 68% confidence region, while the outer lighter-shaded contours 
are the 95% confidence region. The horizontal green bands show the peculiar velocity constraints, while the magenta bands 
show the weak lensing constraints. The central blue contours show the constraints possible when the two datasets are analysed 
together. The cross symbol shows the expected parameterization for General Relativity. 



the covariance of correlation functions in 75 patches of 
area 16 dcg^, and estimate the combined covariance in- 
cluding shape noise and large-scale structure covariance 
using the results of [28] . To model the galaxy distribution 
function n{z) for the CFHTLS sample, we follow [24] in 
using the parametric form 



n{z) oc 



+ z^^ 
zB + C ' 



(9) 



where A, B and C are nuisance parameters to be 
marginalised over. We wish to avoid strongly non-linear 
scales in modified gravity, as these cannot be treated 
properly without using N-body simulations; we there- 
fore remove from our analysis the data in the highly 
non-linear regime, namely 9 < 30 arcmins. 



B. Peculiar velocity measurements 

The coherent peculiar motions of galaxies can be sta- 
tistically estimated through the measurement, and mod- 
eling, of the large-scale redshift-space distortions [29]. 
The radial peculiar velocities can be expressed as a one- 
dimensional velocity dispersion, discussed and de- 
fined in [22, 30], for which 



^ J Pee{k,a)dk , 



(10) 



where Pee is the three-dimensional power spectrum of 
= 6/aH, and has units of h~^Mpc. In [22] the 
one-dimensional velocity dispersion of galaxies was de- 
termined from the two-dimensional two-point correla- 
tion function of SDSS DR7 LRGs, and was found to have 



the values (t„ = B.Ol^o'le ^ ^ Mpc at a mean redshift of 
z = 0.25, and = 3.69lg'47 Mpc * mean red- 
shift of z = 0.38. We will therefore use these values of 
sigmtty as our PV observable. We should note that these 
measurements assume a flat ACDM background, i.e. a 
dark energy equation of state fixed at w = —1. This is 
a valid constraint in our study, since the background ex- 
pansion is known to be close to ACDM, and generically 
a modified gravity can mimic this expansion while hav- 
ing a different structure growth history. Therefore the 
constraints on modified gravity in this study come from 
structure growth. 



C. Current constraints on /i^ and 

In this section, we present the current constraints on 
the MG parameters discussed in Section II from the avail- 
able WL and PV data discussed above. In detail, we vary 
the following set of parameters. 



(11) 



where Wb = flbh^ and cUc = flch^ arc the physical baryon 
and cold dark matter densities relative to the critical den- 
sity respectively, 9^ is the ratio (multiplied by 100) of 
the sound horizon to the angular diameter distance at 
decoupling, r denotes the optical depth to re-ionization. 
Us and As are the primordial power spectral index and 
amplitude respectively, and Eg and are the MG pa- 
rameters for our scale-independent parametrization from 
Eqn. (6). We also vary, and marginalize over, several 
nuisance parameters denoted by Af when performing our 
likelihood analysis, including three parameters (A, B, C) 
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s=l 


s=3 


Ml 


Si 




S3 


WL 


[-6.3,1.9] 


> -0.4 


Unconstrained 


0.32 ± OASto^gi 


PV 


0.06 ±0.20 ±0.40 


Unconstrained 


-0.03 ± O.SO^i;;? 


Unconstrained 


PV+WL 


-0.002 ±0.20 ±0.39 


0.11 ±0.19--- 


-O.OSto TO ± 1-5 


-0.28-0.54 



TABLE I: We present the mean, 68%, and 95% confidence limits for the modified gravity parameters shown in Eqn. 6 for both 
s = 1 and s — 3. We provide constraints based on using the weak lensing (WL) and pecuUar velocity (PV) data alone, as well 
as the combination of these two measurements (PV±WL). We only provide 95% confidence limits for the WL-only constraints 
on fii and Ei. 



for WL data, associated with the galaxy distribution n{z) 
in Eqn. (9), and one for supernovae, which accounts for 
the cahbration uncertainty in measuring the supernova 
intrinsic luminosity. Finally, we assume a flat ACDM 
background cosmology (as discussed in section B). 

Given the set of cosmological parameters P in 
Eqn. (11), wc calculate the expected observables includ- 
ing the CMB shift parameters [31], the luminosity dis- 
tance for supernovae, the growth factor for PV (0) and 
the E-mode component of the weak lensing shear {£,e) 
using MGCAMB [6]. We then constrain the model param- 
eters by comparing these predictions with the PV and 
WL data discussed above, as well as available SNe from 
the UNION-2 sample [32] and CMB shift parameters de- 
rived from WMAP seven year data [31], using a ver- 
sion of the Markov Chain Monte Carlo (MCMC) package 
CosmoMC [33, 34] modified to include our extra parame- 
ters. 

The left panel of Fig. 1 shows our constraints on Ei 
and III for the s = 1 model. In this particular model, any 
departures from the GR value oi fj, = would appear at 
relatively high redshift and hence the effect on the density 
contrast 6 is more prominent than for the s = 3 model. 
One can see a strong degeneracy between fii and Si in 
the WL constraints: if ^1 is large, the growth rate is 
enhanced, which can be compensated for by decreasing 
El. However, this degeneracy does not continue if Ei 
drops below —1; since the WL observables depend only 
on the square of E, a large negative value of Ei will 
actually enhance the lensing signal and can no longer 
compensate for enhanced growth due to large fxi. For 
negative ^1, the growth rate is suppressed and this can 
be compensated for by increasing Ei. In this case there 
is no upper bound for Ei, and no lower bound for fii. 

The right panel of Fig. 1 shows the constraints on the 
parameters for the s = 3 model. In this model, fi is 
now modified at relatively low redshifts, so the effect of 
varying fi^ on the lensing measurements is weaker than 
in the s — 1 model above. As a consequence, the WL 
constraints now become more vertical in the E3— /i3 plane 
as essentially only E3 determines the WL signal. 

Overall, Fig. 1 shows the usefulness of using PV to 
break the E — ^ degeneracy inherent in using just the WL 
measurements. For both models considered, the com- 
bined (blue) PV+WL contours are substantially tighter 
than those obtained from using either the PV or WL 
measurements alone. 



In Table I, we present our measurements of Es and 
fis for the combined WL and PV analysis, as well as 
measurements using the datasets separately. Again, we 
see significant improvement in the constraints on E^ and 
fis when the two datasets are used together as illustrated 
in Fig. 1. For example, using both PV and WL data 
together, we obtain /ii = -0.002 ± 0.59, Ei = O.lllJ] j| 

(for s = 1) and fj,3 = -O-OStl fa^ ^3 = 0.17tJ^ (for s = 
3) at the 95% confidence level, which is fully consistent 
with the expectation for GR (i.e., ^1 = /X3 = Ei = E3 = 
0) although the errors are large, especially for the cubic 
model. 



IV. FUTURE FORECASTS 

In this section, we provide joint error forecasts on Eg 
and fis from a future space-based DE mission such as 
Euclid [35]. 



A. Fisher matrix analysis 

The shear two-point functions and PV measurements 
depend not only on the MG parameters (E^,/!,,) but 
also on the background cosmological parameters and the 
primordial power spectrum. We do not assume that 
these quantities are known, but instead assume that 
CMB data and supernova Type la (SNIa) measurements 
will be available to help constrain them. For the CMB 
power spectra, we include in our analysis the (unlensed) 
Ci spectra of temperature-temperature, temperature- 
polarization and polarization-polarization and use the 
expected errors for the Planck survey [36]. We include 
constraints from a deep space-based SNIa experiment 
[37] with observations of 3000 supernovae with median 
z ~ 0.8, and assume that the SNIa evolution uncertainty 
is well understood. In the forecast, we treat (wf,, lOc, 6s, 
Zreiom "-sj ^s-Wq, Wa, Hsj ^s) ^s frcc parameters with 
mean values of the WMAP7 best fit where the dark en- 
ergy equation of state is modeled as w{a) = wo+Wa{l—a) 
and Zreion is the re-ionization epoch. 

To calculate the expected errors on these parameters, 
we make a first-order Taylor expansion for the param- 
eter dependences of all the observables (CMB, SN, cos- 
mic shear two-point functions and PV). In this "linear 
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FIG. 2: Contour plots of Ss and jis calculated from our Fisher forecast analysis. The left panel shows the s = 1 case, while 
the right panel shows the s — 3 case. Solid curves represent the scenario in which all other cosmological parameters except Es 
and fis are fixed, while the dashed curves represent the scenario in which w = — 1 is assumed, and the dotted curves represent 
the scenario where all the cosmological parameters are allowed to vary including w. 



response" approximation, given the expected experimen- 
tal errors on the povifer spectra, we can easily calculate 
the expected error covariance matrix as the inverse of 
the Fisher matrix. This "linear response" approximation 
(using the first order Taylor expansion) can be improved 
with a careful choice of the parameters including Qs [38]. 

We consider a future WL survey such as Euclid [35] 
covering 20,000 square degrees with sufficient optical and 
infra-red sensitivity to yield a galaxy density of 35 galax- 
ies per square arcminute [39]. We also assume five red- 
shift bins, with a spacing of Az = 0.4, and an underly- 
ing redshift distribution of dn/dz (x z^ exp [— (2/0.64)^'^] 
[39] . The contribution to a Fisher matrix from the shear- 
shear correlations is then given by 



E ■ 

l,a,b,c.d 



21 



(12) 



where C^''' is the angular cross-power spectrum between 
redshift bin a and b, W is the inverse of the total conva- 
riance matrix, and the subscript p denotes differentiation 
with respect to a particular cosmological parameter [40] . 



Using spectroscopic redshift measurements, we are able 
to isolate the PV power spectra from the redshift-space 
power spectrum P°'^''(k) of a galaxy redshift survey. The 
latter is commonly modeled as [41] 

P;'^^(k) = [P^^(k) + 2c2p^e(k) + c^Pee(k)] 

X F{k'c^al{z)), (13) 

where 8 = 9/aH, c is the cosine of the angle to the line of 
sight, Pgg is the true galaxy power spectrum. Pee is the 



peculiar velocity power spectrum and Pge is the galaxy- 
velocity cross spectrum. The separation of Pgg(k.) and 
Pee(k) is made possible using the angular dependence 



of P°''''(k). 



In order to obtain cosmological constraints from the 
peculiar velocities, we need to understand the error co- 
variance matrix of the observables. We estimate this us- 
ing the Fisher matrix [42] 



Tiobs 



{ki, Zj) 



k^dk 
2(2^ 



dc Voff(fc,c, Zj) 



d In Pf"" (fc, c, Zj ) d In P°'°''{k, c, zj ) 



dpa 



(14) 



where a and /? both run from 1 to 2 and denote Pgg and 
Pee respectively. Note that in the linear regime P^e is 
simply the square root of the product of Pgg and Pee^ 
so we take this to be the case in this work. The effective 
volume V^ff in each redshift bin j is given by 



Ve({{h,c,Zj) 



niP°' 



nJP°bB(fc^,c,Zj) + l 



survey 



{z= 



(15) 

where is the shot noise term coming from the fi- 
nite galaxy density, and Vsmvcy {zj) is the survey vol- 
ume in a given redshift bin. On large scales, the cosmic 
variance term dominates over the shot noise term and 
Vcs{ki, ^, Zj) is nearly identical to Vsurvcy{zj)- For our 
estimation, we consider a 20000 sq deg Euclid-like sur- 
vey with H-a emitter number density as a function of 
redshift given by [43] (Table 2, limiting flux 4 x 10"^^ 
erg s~^ cm~^). 

Following the above estimation of Pgg , Pee ^^nd their 
errors, we can then constrain the cosmological parame- 
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s = 1 
cr(Ei) (7{fii) 


s = 3 
o-(E3) (7(^3) 


Current constraint 


0.59 0.59 


0.88 2.3 


Future forecast 


Fix cosmological params. 
Assume only k; = — 1 
Allow w to vary 


0.0081 0.0046 
0.035 0.041 
0.037 0.046 


0.017 0.014 
0.040 0.069 
0.066 0.12 



TABLE II: The predicted errors on Es and Hs for a future Euclid-like satellite mission compared to the current constraints 
obtained in this paper. 



ters from the decomposed Pee ■ (Here we elect not to use 
Pgg on account of the bias). The Fisher matrix for cosmo- 
logical parameters, using the peculiar velocity spectrum, 
can be written as [44] 



pp 



E 



dPee 



ee 



dp (P-l°b^)22(fez,Zj) dp' 



(16) 



where i denotes k bins up to fc = 0.1/iMpc ^, and j 
denotes redshift bins from z = to 2 with spacing Az — 
0.2. 



B. Future forecasts 

In this subsection, we study the constraints on the MG 
parameters in three different scenarios: 1) we only vary 
the parameters S^, ^s, fixing all other cosmological pa- 
rameters; 2) we vary all cosmological parameters with 
the assumption of a fiat ACDM expansion history; 3) 
we vary all cosmological parameters including w using 
w{a) = wq + Wa{l — a). In each case, we again consider 
the s = 1 and s = 3 models. 

Fig. 2 shows constraints on and fig in the three 
scenarios above. Again, thanks to the complementar- 
ity of WL and PV measurements, the constraints on 
these MG parameters are very tight (solid contours) if 
we fix the other cosmological parameters. However, if we 
marginalise over the uncertainties in the other cosmolog- 
ical parameters, then the constraints on Es and fis are 
degraded by an order of magnitude: the dotted contours 
in Fig. 2 represent the full marginalization over other cos- 
mological parameters allowing w to vary, while the inner 
dashed contours represent the w = —1 case. 

In Table 2, we present the expected constraints on 
and for the three scenarios discussed above. For com- 
parison, we also provide in the table the current con- 
straints on these MG parameters and even in the most 
conservative case (allowing w to vary), the future satel- 
lite constraints should be at least 20 times better. It is 
interesting to note the difference between the 5 = 1 and 
the s = 3 model. As s decreases, the expected departure 
from GR starts at an earlier epoch thus resulting in a 
stronger constraint. 

Comparing the current constraints in Fig. 1 and the 
future forecasts in Fig. 2, we notice a change in the ori- 
entation of the Es — /is contours. Namely, in the current 
constraints, there is an anti-correlation between S. and 



fis, while for the future forecasts, this turns into a pos- 
itive correlation. We note that without inclusion of the 
PV information there is still an anti-correlation between 
Es and /is in the forecast as well, so it is the inclusion of 
PV which is changing the orientation. 

This change in orientation is due to the effect of 
marginalising over other cosmological parameters. To 
better understand this effect, let us examine the effect 
of marginalizing over cOm- When we examine the covari- 
ance between different cosmological parameters for a PV 
experiment, we find that both Es and /is are strongly 
anti-correlated with LUm, and to a lesser extent with wq, 
Us, Os, LiJb and As] this means, for instance, that an in- 
crease in /i can be compensated for by a decrease in ujm- 
Equally, if one decreases Wm, the WL data can be fit- 
ted if one also increases E. As a result, given precise 
measurements for both PV and WL, Eg and /is can be- 
come positively correlated, as they may both change in 
the same sense to oppose a change in lo„i in the oppo- 
site sense. However, if the PV data are sufficiently noisy, 
the anti-correlation of /i and E exhibited in the WL con- 
straints alone (see Fig. 1) will dominate. 



V. TOWARDS MODEL INDEPENDENT 
CONSTRAINTS ON MG PARAMETERS 

As we saw in section HI and section IV, constraints 
on /is (and to a lesser degree Es) depend strongly on the 
assumed value of s. One may wonder if it is possible 
to extract any model-independent information about /i 
from the current PV and WL data. In principle, one can 
bin /i into a large number of narrow redshift bins and 
find their best constrained uncorrelated linear combina- 
tions - the so-called eigenmodes [7, 45, 46]. In practice, 
this would be a numerically challenging project, given the 
large range of z over which variations of /i can affect the 
growth. However, one would expect that with PV and 
WL measurements at a single redshift z*, one would con- 
strain just one of the eigenmodes, corresponding to a net 
change in /i over Q < z < z* . The tightness of the con- 
straint would depend on how early the departures from 
/i = 1 are allowed to start, because an earlier change in /i 
results in a larger change in the growth factor throughout 
z. 

To demonstrate this point, we try fitting /i (which can 
be any constant over a wide bin) between z^ and z = 0, 
and then setting /i = 1 for z > z^. Note that this is not 
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quite the same as the s = case in Eqn. (6), as in this 
equation we do not cxphcitly require fj, to be exactly unity 
above a certain redshift. The results for the choices of 
Zs = 2 and = 10 are shown as shaded boxes in the top 
right and left panels respectively of Fig. 3. As expected, 
the constraint on the constant bin is much stronger in 
the Zs = 10 case. 

However, the analysis of the SDSS DR7 LRGs pro- 
vides measurements of the peculiar velocities of galaxies 
at more than one redshift, i.e. we have measurements at 
zi ~ 0.25 and Z2 = 0.38. Therefore, one could ask if it is 
possible to obtain additional information about the vari- 
ation of fi between Z2 and zi independent of the assumed 
value of Zs- To address this, we examine ^ in two bins, 
namely fiA in [zi,Z2] and fis in [^2,2:3], with /i = 1 for 

Z > Zs. 

The constraints on these two bins are shown again in 
Fig. 3 for the two cases of Zs = 2 and Zs = 10. As 
expected, the constraint on fi is tight in the higher red- 
shift bin, with the error bar being much smaller for the 
Zs = 10 case. On the other hand, the errors on the low 
redshift bin are much larger. Interestingly, however, they 
are effectively independent of the choice of Zs , as is the 
correlation coefBcient between the two bins. 

By diagonalising the covariance matrix we find the un- 
correlated linear combinations of the two bins to be = 
0.9973^A - 0.0737/iB and qb = 0.0737yU^ + 0.9973yUB for 
the Zs = 2 case and qa ~ 0.9993^a — 0.0364/iB and 
qb = 0.0364yU^ + 0.9993yU_B for the Zg = 10 case. One can 
see from Fig. 3 that the error on the better constrained 
eigenmode depends strongly on Zs, while the error on 
the second eigenmode stays the same. This demonstrates 
that having measurements at multiple redshifts can pro- 
duce constraints on the variation of fi at low redshifts 
that are independent of its assumed time-dependence at 
high redshifts. Finally, we note that the values for /i 
found in the two bins do not show any deviation from 
the GR predictions. 



VI. CONCLUSION 

In this paper, we have explored the complementarity 
of weak lensing and peculiar velocity measurements to 
obtain better constraints on the modified gravity param- 
eters fi and E, as defined in Eqns. (4) and (5). Peculiar 
velocity measurements are only sensitive to variations of 
/I, thus breaking the degeneracy between /i and S found 
when using just the weak lensing measurements. Using 
a simplified parametrization for ^ and E (see Eqn. 6), 
we have derived the first constraints on fj, and S using 
weak lensing measurements and peculiar velocity mea- 
surements (from CFHTLS and SDSS respectively). As- 
suming a flat ACDM background cosmology, so that here 
any modifications to gravity are probed by structure 
growth, the best constraints we find using both datasets 
together are m = -0.002 ± 0.59 and Ei = O.lltg;'^^ 
for s = 1, and /13 = -0.08^2:28 and E3 = OAltoH 











H 


1- 






' 

h-i^ - 




% 


% 


H 
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q - 0.9993^1^- 0.0364^^^ 
q- 0.0364H + 0.9993|i„ 
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FIG. 3: Constraints from PV data on as a function of 
redshift. We show the results of using 2 redshift bins, with 
fi = fiA within [21,22] and /i = /is witliin [22,2s]. It is as- 
sumed that mu = 1 above 2s = 2 (left) and above Za = 10 
(right). In addition we display (grey box) results for constant 
fi below redshift 2s = 2 (left) and below Zs = 10 (right) when 
PV is combined with CMB shift parameter and SNe. The 
lower panels show the linear combination of ha and fiB that 
diagonalise the covariance matrix. 



for s = 3 (95% confidence limits). These results are 
consistent with the expectation from General Relativity, 
namely fii = /i3 = Ei = E3 = 0. As shown in Table 1, 
these constraints arc much worse if the WL and PV data 
are used separately, in which case some of the parameters 
become unconstrained. 

We also performed a Fisher error forecast for a space- 
based Dark Energy mission like Euclid. For example, 
assuming a flat ACDM cosmological background, we pre- 
dict that the constraints on the modified gravity pa- 
rameters are improved by at least an order of magni- 
tude or more, e.g. cr(/ii) = 0.037, o'(Ei) = 0.032 and 
aifis) = 0.068, cr(E3) = 0.038 (see Fig. 2 and Table II 
for the full details). It is interesting to note that the 
expected precision on these errors is comparable in size 
to the present-day uncertainties on cosmological param- 
eters like w (approximately known to 10% today) and 
(~ 5% error from present observations). Therefore, 
a Euclid-like mission will deliver outstanding constraints 
on modified gravity as well as improving errors on the 
background cosmological parameters. 

That said, we found that the constraints on fi depend 
strongly on the assumed time variation of this parameter, 
i.e. s in Eqn. (6). We have therefore proposed a method 
for extracting model-independent information about fi 
by examining fj. in different redshift bins. For the SDSS 
DR7 PV measurements, we have used two bins and are 
able to obtain a model-independent constraint on fi at 
low redshifts which is again consistent with expectations 
for GR (see Fig. 3). This demonstrates that having PV 
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measurements at multiple redshift intervals can produce 
a constraint on the variation of /x, at low redshifts, which 
is independent of the assumed time dependence at high 
redshifts. This technique will become increasingly im- 
portant in the next few years as new peculiar velocity 
measurements are published from surveys like WiggleZ 
[47], BOSS and VIPERS. 
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